Wigner-function theory and decoherence of the quantum-injected optical parametric 

amplifier 
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Recent experimental results demonstrated the generation of a quantum superpositon (MQS), 
involving a number of photons in excess of 5 x 10 4 , which showed a high resilience to losses. In order 
to perform a complete analysis on the effects of de-coherence on this multiphoton fields, obtained 
through the Quantum Injected Optical Parametric Amplifier (QIOPA), we invesigate theoretically 
the evolution of the Wigner functions associated to these states in lossy conditions. Recognizing 
the presence of negative regions in the W-representation as an evidence of non-classicality, we focus 
our analysis on this feature. A close comparison with the MQS based on coherent \a) states allows 
to identify differences and analogies. 

PACS numbers: 



I. INTRODUCTION 



In the last decades the physical implementation of 
macroscopic quantum superpositions (MQS) involving a 
large number of particles has attracted a great deal of 
attention. Indeed it was generally understood that the 
experimental realization of a MQS is very difficult and in 
several instances practically impossible owing to the ex- 
tremely short persistence of quantum coherence, i.e., of 
the extremely rapid decoherence due to the entanglement 
established between the macroscopic system and the en- 
vironment P, 0, [H, HI- Formally, the irreversible decay 
towards a probabilistic classical mixture is implied theo- 
retically by the tracing operation of the overall MQS state 
over the environmental variables [H, Q . In the framework 
of quantum information different schemes based on op- 
tical systems have been undertaken to generate and to 
detect the MQS condition. A Cavity - QED scheme 
based on the interaction between Rydberg atoms and 
a high-Q cavity has lead to the indirect observation of 
Schrodinger Cat states and of their temporal evolutions. 
In this case the microwave MQS field stored in the cavity 
can be addressed indirectly by injecting in the cavity, in 
a controlled way, resonant or non-resonant atoms as ad 
hoc " measurement mouses" @, H| ■ A different approach 
able to generate freely propagating beams adopts photon- 
subtracted squeezed states; experimental implementa- 
tions of quantum states with an average number of pho- 
tons of around four have been re por ted both in the pulsed 
and continuous wave regimes [gL llOL [Til . [HJ . These states 
exhibit non-gaussian characteristics and open new per- 
spectives for quantum computing based on continuous- 
variable (CV) systems, entanglement distillation proto- 
cols [HI [3] , loophole free tests of Bell's inequality. 

In the last few years a novel "quantum injected" opti- 
cal parametric amplification (QI-OPA) process has been 
realized in order to establish the entanglement between 
a single photon and a multiphoton state given by an av- 
erage of many thousands of photons, a Schrodinger Cat 



involving a " macroscopic field" . Precisely, in a high-gain 
QI-OPA " phase-covariant" cloning machine the multi- 
photon fields were generated by an optical amplifier sys- 
tem bearing a high nonlinear (NL) gain g and seeded 
by a single-photon belonging to an EPR entangled pair 

[StmBHa-tiij. 

While a first theoretical insight on the dynamical fea- 
tures of the QIOPA macrostates and a thorough experi- 
mental characterization of the quantum correlations were 
recently reported (20l . [2l| . a complete quantum phase- 
space analysis able to recognize the persistence of the 
QI-OPA properties in a decohering environment is still 
lacking [22|, |23|. Among the different representation of 
quantum states in the continuous- variables space [24| , the 
Wigner quasi-probability representation has been widely 
exploited as an evidence of non-classical prop erties, such 
as squeezing [25[ and EPR non-locality |26| . In partic- 
ular, the presence of negative quasi-probability regions 
has been considered as a consequence of the quantum 
superposition of distinct physical states (27| . 

In the present paper we investigate the Wigner func- 
tions associated to multi-photon states generated by op- 
tical parametric amplification of microscopic single pho- 
ton states. We focus our interest on the effects of de- 
coherence on the macro-states and on the emergence of 
the "classical" regime in the amplification of initially pure 
quantum states. The Wigner functions of these QI-OPA 
generated states in presence of losses are analyzed in com- 
parison with the paradigmatic example of the superposi- 
tion of coherent, Glauber's states, |a). 

The paper is structured as follows. In Section II, we in- 
troduce the conceptual scheme and describe the evolution 
of the system both in the Heisenberg and Schrodinger 
picture. Section III is devoted to the calculation of the 
Wigner function of the QI-OPA amplified field. We first 
consider a single-mode amplifier, which is analogous to 
the case of photon-subtracted squeezed vacuum. Then 
we derive a compact expression of the Wigner function 
in the case of a two- mode amplifier in the " collinear" case, 
i.e. for common k— vectors of the amplified output fields. 
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In Section IV, we introduce, for the collinear de- 
coherence model apt to simulate the de-cohering losses 
affecting the evolution of the macrostates density matrix. 
This evolution is then compared to the case of the coher- 
ent | a) MQS. Section V is devoted to a brief review of the 
features of coherent states superpositions (CSS). Hence 
in Section VI we derive an explicit analytic expression of 
the Wigner functions for the QI-OPA amplified states in 
presence of decoherence. The negative value of the W- 
representation in the origin of the phase-space for such 
states above a certain value of the " system-enviroment " 
interaction parameter, is an evidence of the persistancc 
of quantum properties in presence of decoherence. This 
aspect is then compared with the case of the \a) states 
superposition. Then in Section VII we study a comple- 
mentary approach to enlighten the resilience to losses of 
the QI-OPA MQS with respect to the \a) states superpo- 



sition. Precisely, we define a coherence parameter, based 
on the concept of state-distance in Hilbert spaces, and we 
study its decrease as a function of the losses parameter for 
both classes of MQS. Finally, in Section VIII we analyze 
a different configuration based on a non-collinear opti- 
cal parametric amplifier in the quantum injected regime, 
i.e. a universal quantum cloning machine. We calculate 
the Wigner function associated to the states generated 
by this device in absence and in presence of decoherence, 
focusing on the persistence of quantum properties after 
the propagation over a lossy channel. 



II. OPTICAL PARAMETRIC AMPLIFICATION 
OF A SINGLE PHOTON STATE 




FIG. 1: Scheme of the experimental setup. The main UV laser beam provides the OPA excitation field beam at A = 397. 5nm. 
A type II BBO (Beta Barium Borate) crystal (crystal 1: CI) generates pair of photons with A = 795nm. In virtue of the 
EPR non-local correlations established between the modes Ica and ks, the preparation of a single-photon on mode ks with 
polarization state is conditionally determined by detecting a single-photon after proper polarization analysis on the mode 
kA (polarizing beamsplitter (PBS), A/2 and A/4 waveplates, Soleil-Babinet compensator (B), interferential filter (IF)). The 
photon belonging to ks, together with the pump laser beam h' p , is fed into an high gain optical parametric amplifier consisting 
of a NL crystal 2 (C2), cut for collinear type-II phase matching. The fields are coupled to single-mode fibers. For more 
details refer to .21]. Right Inset: Measurement apparatus. After fiber polarization compensation, the field is analyzed by 
two photomultipliers {P+,P_} and then discriminated through an O-Filter device. Left inset: Action of the O-Filter in 
the photon number space: the (±1) outcomes are assigned whether — m 7r± > k or m 7r± — > k, where k is a tunable 
threshold condition. The central region leads to an inconclusive (0) outcome, and the two orthogonal macrostates cannot be 
discriminated. 

I 



As a first step we consider the generation of a mul- tiphoton quantum field, obtained by parametric am- 
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plification. Let us briefly describe the conceptual 
scheme. An entangled pair of two photons in the 



{\H) A \V) B -\V) A \H) B ) is 
down- 



singlet state \ip iab 

produced through a spontaneous parametric 
conversion (SPDC) by crystal 1 pumped by a pulsed UV 
pump beam: FigfT] There \H) and |V) stands, respec- 
tively, for a single photon with horizontal and vertical 
polarization i/rt) while the labels A, B refer to particles 
associated respectively with the spatial modes k^and kp. 
The photon belonging to kg , together with a strong UV 
pump beam, is fed into an optical parametric amplifier 
consisting of a non-linear crystal 2 pumped by the beam 
kp. The crystal 2 is oriented for "collinear operation", 
i.e., emitting pairs of amplified photons over the same 
spatial mode which supports two orthogonal ~tt modes, 
respectively horizontal and vertical. Let us analyze the 
properties of the resultant amplified field. 



A. Collinear optical parametric amplifier 



The complete Hamiltonian of the system reads: 



H 



Huj (a^an + a\,av^ + i^X (^li^v e 



ancive 



where x is proportional to the 2 nd order non-linear sus- 
ceptibility and to the amplitude of the pump field Eq ; 
a H and a,y are the creation operators associated to the 
mode k respectively, with polarization w horizontal and 
vertical, H, V. We assume a classical and undepleted co- 
herent pump field. Let us consider only the interaction 
contribution to H and a reference system which rotates 
with angular speed 2u> . The time-independent interac- 
tion Hamiltonian is found to be: 



Hi = ihx — anctv 

The Heisenberg evolution equations are: 



da v 1 
dt ih 

with general solution: 



a n , Hi 



(2) 



(3) 



a n (t) = a n (0) cosh (xt) + a[ ± (0) sinh (xt) (4) 

where tt = H, V. 

The average number of photons generated in the SPDC 
process can be easily calculated in the Heisenberg picture 
formalism by applying the evoluted operators to the ini- 
tial vacuum state \0H,0V) = |0). Heretofore, \nir,rrnr±) 
stands for the Fock-state with n photons with tt polar- 
ization and m photons with the orthogonal tt± one. For 
the horizontal polarization we get: 



(0(4, (t) a H (t) |0) = sinh 2 ( X t) = m 
The same result holds for the vertical polarization. 



(5) 



B. Output wavefunctions 

In any "equatorial" polarization basis {\ip) , j^r 1- ^}, re- 
ferred to a Poincare sphere with "poles" \H) and \V), 
the Hamiltonian of the polarization non-degenerate opti- 
cal parametric amplifier, described by the expression ((2|), 
can be expressed as: 



Hi = 



((a*) s 



<p±J 



(6) 



where the corresponding field operators are: 



(a% + e^a f y ) , fit ± = _i_ ^ _ c . ^ 

The Hamiltonian can then be separated in the two po- 
larization components {7? v ,7? v i}. The time-dependent 
field operators are: 

a v (t) = 6 V (0) cosh(xt) + e-^al(0) sinh( x i) (7) 

Let us now restrict the analysis to the basis {7? + ,7f_}, 
in which we calculate the output wavefunction in the 
Schroedinger picture; the unitary evolution operator is 



U = exp 



with 



.Hit 



C/ ± =exp[±!(V ± 2 -a 2 ± )] 



U+U- 



(8) 



(9) 



where g = xt is the non linear gain of the process. 

The expression of U enlightens the decoupling between 
the two polarization modes. A simple expression of the 
operators U+ and £/_ can be obtained adopting the fol- 
lowing operatorial relation [28| : 



U± 



exp 



±r 



s/2 



exp 



-ln(cosh(7) ( a^ ± a± + — 



x exp 



a± 
V2 



(10) 



with r = tanh g. Let us consider the injection of a single 
photon state with generic equatorial polarization Tr v into 
the OPA: 

IV>i«> = \<p) = -±=(\H) + e i *\vj) = 

= e% (cos | |l+,0-} +isin| |0+, 1-}) (11) 

The multiphoton output state of the amplifier is found: 



1 4-, 



I <!>'"> ^ ( OS | |$+) 



/hin-|$~) (12) 



where 
s± 



* ± ) = 1>±L> t |1±,0t> = (l>±|i> ± ) (u t \o) t ) (13) 



4 



Simple calculations lead to the expressions: 



CP 



E 



/2JT v/(2i + 1)! 



(2i + l)±,2j T ) 
(14) 



where C = cosh 5. The quantum states |$ + ) and |<&~) 
are orthogonal, being the unitary evolution of initial or- 
thogonal states. We observe that |<£> + ) presents an odd 
number of 7? + polarized photons and an even number of 
7r_ polarized ones. Conversely for |$~). 

The average number of photons with polarization 7? + 
can be estimated in the Heisenberg representation as: 

(h + ) = (ip in \ a\_(t)a+(t) \ip in ) — fh -f~ (2m + 1) cos 2 {^\ 

(15) 

where m = sinh (g). The phase dependence shows that 
for ip = 7r the average number of generated photons is 
equal to the spontaneous one, while for ip = 0, corre- 
sponding to the stimulated case, an increase of a factor 3 
is observed in the average number of 7? + polarized pho- 
tons. 

Likewise, the average number of photons 7?_ polarized 
is given by: 



(h-) = (ip in \ al(t)a-(t) \ip in ) = m + (2m + 1) sin 2 (|) 

Hence, by varying the phase ip we observe a fringe pattern 
which exhibits a gain dependent visibility: 



X i) _ (n+ (<p = 0)) - (n- (p = 0)) 2m + 1 

" l (»+ (ip = 0)) + <n_ {<p = 0)) 4m +1 1 ' 



v th 



In the asymptotic limit (g — > 00) 

We now briefly discuss the phase- covariance proper- 
ties of the optical parametric amplifier when injected 
by a single-photon state. In this configuration, this de- 
vice acts as optimal phase- covariant cloning machines, 
and due to the unitary evolution of the process, the su- 
perposition character of any generic input state \<p) = 
^= (\H) + e lv \V}) is maintained after the amplification: 



V2 



o*<P\3> v 



(18) 



Let us stress that there is freedom in the choice of the 
macrostates basis vectors. For example, the l^*) states 
can be used in the expansion of the overall state as in 
(fT2"| . Hence, the equatorial |<!> v ) amplified state can be 
written as the macroscopic quantum superposition cither 
of the {|$+),|$-)} and the {\$ H ), "basis" states. 

Furthermore, due to their resilience to decoherence [29| . 
all equatorial macro-qubits represent a preferred " pointer 
states basis" Q for writing in the form of a macro- 
scopic quantum superposition. 



III. WIGNER FUNCTIONS OF THE 
AMPLIFIED FIELD 

In order to investigate the properties of the output field 
of the QI-OPA device in more details, we analyze the 
quasi-probability distribution introduced by Wigner [23| 
for the amplified field. The Wigner function is defined as 
the Fourier transform of the symmetrically- ordered char- 
acteristic function x(v) of the state described by the gen- 
eral density matrix p 



X(V) = Tr [pexp (jia) - if a)] 



(19) 



The associated Wigner function 

W(a) = i J exp [rj*a -r)a*)x (v) A (20) 

exists for any p but is not always positive definite and, 
consequently, can not be considered as a genuine prob- 
ability distribution. Since the early decades of quantum 
mechanics, the emergence of negative probabilities has 
been identified as a peculiar feature of quantum physics, 
and has been connected to the mathematical segrega- 
tion of states which physically live only in combination 
[27l ]. In parallel, the non classicality of a quantum state 
is expressed by a Glauber's P-representation [Hi Ull HH 
which is more singular than a delta function, i.e. the 
P{a) proper of coherent states. This means that the 
system does not possess an expansion, in terms of the 
overcomplete semi-classical \a) state basis, that can be 
interpreted as a probability distribution. However, the 
negativity of the Wigner function is not the only param- 
eter that allows to estimate the non-classicality of a cer- 
tain state. For example, the squeezed vacuum state [25[ 
presents a positive IF-representation, while its proper- 
ties cannot be described by the laws of classical physics. 
Furthermore, recent papers have shown that the Wigner 
function of an EPR state provides direct evidence of 
its non-local character [2(1 [32j , while being completely- 
positive in all the phase-space. 

The complex variable a in (J2TJJ) is the eigenvalue asso- 
ciated to the non-Hcrmitian operator a which acts upon 
the coherent state \a) as follows: a\a) = a\a). It is pos- 
sible to decompose a in its real and imaginary part and 
then to define the quadratures operators X and Y which 
allow the representation of the field in the phase space 
a = X ~*2 ■ The quadratures operators are Hermitian 
operators and thus correspond to physical observables 
proportional to the position q and the momentum p fol- 
lowing the relations: X = \f^-q and Y = —\J~^P- The 
uncertainty principle leads to: AX AY > 1. 



A. Single mode amplifier 

For the sake of simplicity, let us first consider the 
Hamiltonian of a degenerate amplifier acting on a sin- 
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gle k— mode with polarization 7T+: 
Hi 



in— (ij 



(21) 



When no seed is injected, the amplifier operates in the 
regime of spontaneous emission and the characteristic 
function reads: 



As a further example, we consider the injection of the 
2-photon state \ipi n ) = |2+). We obtain 

X2(V, t) = (l - 2 Ht)\ 2 + \ |r,(t)| 4 ) exp \g(t)^ (28) 



Xofat) = (0|e»pfo4(t) -»?*&+(*)] |0) 

= (o| exp^wa^ - *i*(t)a+] jo> 



and the Wigner function reads: 
2 

TV 



(22) 



with a+(t) = a+ cosh(g) + a\_ sinh(g), r](t) = ?7Cosh((7) — 
rf sinh(g), and g = yt. Hereafter, we explicitly report the 
dependence of the Wigner function from the interaction 
time t. We obtain, using the operatorial relation exp(A + 
B) = expiexpBexp(-l/2[i,i?]): 



Xo fat) = exp 



-l\v(t)f 



(23) 



The calculation then proceeds as follows. Starting from 
the definition of the Wigner function (|2T)|) , we perform the 
two subsequent transformations of the integration vari- 
ables cPr) —> d 2 rj(t) — > xdxd(p, where n — ► rj(t) has been 
defined previously and can be expressed as rj(t) = xe %v . 
The Wigner function then reads: 



W\ 0+) (a,t) 



x[(a—ot*) cos ip-\-i(- 



ol — oi. ) sin 



jI xdxd(p 



W [2+) (X, Y, t) = - [1 - 8 [X 2 e' 2cj + Y 2 e 2g ) 

+2 (X 2 e- 2g + y 2 e 2fl ) 2 ] exp [-2(X 2 e" 2s + yV 9 )] 

As in the single photon case, the negativity of the Wigner 
function is maintained after the amplification process, 
emphasizing the quantum properties of the injected field. 
All these features are shown in Figs. [2] and [3l which 
report respectively the plots of the Wigner functions and 
of their sections X = and Y = for QI-OPA amplified 
states, with the injection of the |0), |1), \ 2) Fock states. 
The analysis of these figures shows that the amplification 
effect leads to the increase of the degree of squeezing for 
the multi-photon output field. The uncertainty on one 
of the two quadratures is decreased while the other one 
is increased coherently with the Heisenberg uncertainty 
principle. Let us note that the quantum character of the 
Fock states is underlined by the negativity of the Wigner 
function in the central region of the quadratures space. 

Finally, this result can be generalized by analogous cal- 
culation to the generic \N+) input state, leading to the 
Wigner function: 



2 f°° ( 1 2 \ 

— / Jo (—2 lal x) exp I ——x I xdx 

n Jo \ 2 / 



■ exp 



-2\a\ 



(24) 

where a = a cosh(g) — a* sinh(g), and Jq(x) is the Bcsscl 
function of order 0. We can now write a = Re(a)e" 9 + 
ilmfaje 9 as a function of the A, Y quadrature operators, 
defined by the expression a = X + lY . By substituting 
such variables Eq. ([24j) becomes 



W {0+) (A, Y,t) = - exp [-2 (A 2 e- 2 « + Y 2 e 2 *)] (25) 

When we consider the case in which a single photon 
with polarization tt + is injected: \ipi n ) = analogous 
calculations leads to the characteristic function 



Xi fat) - (l\exp[r ] (t)al-i 1 *(t)a + }\l) 
l-\v(t)\ 2 ) exp 

The Wigner function reads: 



\m\ 2 



(26) 



2 






1 














7T 





l-4|a| 



exp 



-2I0I 



exp [-2(X 2 e~ 29 + Y 2 e 2g )] (27) 



W ]N+) (a,t) - ~(-l) N L N (4|cf| 2 )e 



-2|«| 2 



(29) 



where Ln are Laguerre polynomials of order N. For all 
N, the non-classical properties of the injected state are 
maintained after the amplification process, as the OPA 
is described by a unitary evolution operator. 

We conclude this Section on the single-mode ampli- 
fier by emphasizing the connection, shown in [33|, be- 
tween the single photon subtracted squeezed vacuum, 
and the squeezed single-photon state. It is found that: 



s 2 



> 2 



aS(O|0) = 5(011), where 5(0 = exp 
is the single-mode degenerate squeezing operator, and 
£ = se is the complex squeezing parameter. This 
evolution operator is obtained by the single-mode am- 
plifier interaction Hamiltonian of Eq. (f2Tj) . For small 
this state possess a high value of overlap [34[ with 
the (\a) — I — a)) quantum superposition. This connec- 
tion between photon-subtracted s quee zed vacuum and 
squeezed Fock-states was extended [33j to the more gen- 
eral p-photon case, obtaining: 



where: 



^ 2 fc fcVb-2fc)! 
x (e 10 sinh s cosh s) k \p — 2k) 



(30) 



(31) 
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WIGNER FUNCTIONS FOR SINGLE MODE OPA 
INJECTED QUBIT AMPLIFIED FIELD 




SECTIONS OF WIGNER FUNCTIONS FOR SINGLE MODE QI-OPA 
INJECTED QUBIT AMPLIFIED FIELD 

x=o 




FIG. 2: Wigner function of the injected states and the corre- 
sponding amplified states, for a single mode degenerate ampli- 
fier with gain g — 4. The axes report the quadratures values 
X and Y. |0) injection of the vacuum state. |1) injection of 
a single photon. |2) injection of two photons. The amplified 
field is plotted on a different scale with respect to the injected 
qubit one, due to the high degree of squeezing introduced by 
the amplification process. 



where Af is an opportune normalization constant. Hence, 
the p-photon subtracted squeezed vacuum is analogous 
to the amplified state of a quantum superposition of 
odd (or even) Fock-states (|3"Tj) . For an increasing value 
of p, the overlap between the c^S 1 ^) |0) states and the 
(|a) ± | — a)) states [35] is progressively higher, but cor- 
responds to the amplification of a more sophisticated su- 
perposition of Fock-states. 



B. Two modes amplifier 

In order to investigate the collinear QI-OPA we have 
to analyze the non-degenerate OPA Hamiltonian ([1]), i.e. 
acting on the both orthogonal polarization modes tth 
and 7?y. For a given input state in the amplifier \tpi n ) 
the characteristic symmetrically- ordered function can be 
written as: 




FIG. 3: Sections of the input states Wigner functions and of 
the relatives amplified fields along the directions X — and 
Y = 0. Since the injected qubit shows a rotational symmetry, 
we report only the X — section. |0) injection of the vacuum 
state. |1) injection of the single photon state. |2) injection of 
the two-photon state. 



where the time dependent operators solve the Heisenberg 
equation of motion ([7]) and are expressed in the basis 

It is useful to rewrite the expression (|32|) by using the 
{7? + ,7?_} polarization basis. We obtain: 



xexp U(t)^-e(t)^=\ |Vin) (33) 



where 

V (t) = (v + 0C-(v* +C)S (34) 
ffl = (V-0C+(V*-C)S (35) 

and C = cosh(<7), S = sinh(g). 

Let us consider the case of single photon injection on 
the polarization mode 7? + , the input state wavefunction 
can be written as \ipi n ) — |l+,0— ). Hence the character- 
istic function is 



XAV,U) = (^ m \e^H(t)-v'a H (t)) e ( S aUt)-e&v(t)) |^ 



XiAv,U) = f 1 - M-j e X p[-i(|r?(t)| 2 + m\ 2 )} 

(36) 

In this case the Wigner function is the quadri- 
(32) dimensional Fourier transform of the characteristic func- 
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tion given by [16| : 



eM-\(Mt)\ 2 + \m\ 2 )] (i-M 



exp -|A| 



2 

7T 

where we have used: 



1 - |A A + Ai 



IAI 



|7A+| 2 + |7A-| 2 + |7B+| 2 + |7B-| 2 



A 



K 



V2 



(37) 

(38) 
(39) 



with K = A,B. 

In this case the squeezing variables are ja+ and 7s_ 
while 7s + and ja- are the respective conjugated vari- 
ables: 



j A + = (a + (3*)e 9 ; -f B + = 
7a- = «(a - /3*)e 9 ; 7s- 



(a* + /?)e~ 9 (40) 
= z(/?-c/)e 9 (41) 



The integration of (|3T|) is analogous to the one reported 
in section UlI AI 

If a — \a\e~xpitp a and [3 = \[3\ exp ipp have a well- 
defined phase relation, ip a = —<pp = ip, then for every ip 
value it is possible to represent the Wigner function in a 
three-dimensional graph, reported in Fig[5J that is, the 
projection of the total function onto a certain subspace. 
The quadrature variables in this subspace are: X = (a + 
/?*) e Y — (J3 — a*). Furthermore, in Fig [5] we report 
the X = and Y = sections of the Wigner function 
for the single photon amplified state in comparison with 
the injected seed. We again note the resilience of the 
negative region, centered in the origin of the phase-space, 
and the presence of the degree of squeezing induced by 
the amplifier. 

When the input state is the state of two photons with 
7r+ polarization \ipi n ) — |2+, 0— ) the characteristic func- 
tion is: 



\v(t)\ 



and the Wigner function is: 



W\ 2+ , -)(a,l3,t) = 



(42) 



exp 



(-|A| 2 ) 



x ( i - 2 |A A + A s | 2 + | |A A + A s | 4 



(43) 



In Fig|4] we report the plots of the Wigner function of 
both the single- and double-photon amplified states com- 
pared to the original W-representations of the injected 
seed. 



INJECTED QUBIT 

|i)Jo) 



WIGNER FUNCTIONS FOR TWO MODES QI-OPA 
AMPLIFIED FIELD 




-0.2 \ 
-0.t\ 




7 0.04 
/ 0.02 



7 -0.02 
-0.04 




100 J _o.04 



FIG. 4: Wigner function of the injected states and the cor- 
responding amplified states, for a non-degenerate amplifier 
with gain g = 4. The axes scales on different graphs are dif- 
ferent due to the high degree of squeezing introduced by the 
amplifier. |1 + , 0—) injection of a single photon 7r+ polarized. 
|2+,0— } injection of a two-photon state tt+ polarized. 



These Wigner functions can also be evaluated by em- 
ploying the results obtained in the single mode amplifier 
case. In this context the characteristic function factorizcs 
into two parts X+iVi t)X— {Wi that refer to polarizations 
7r+ and 7?_ respectively, as a consequence of the indepen- 
dence of the two oscillators. Analogously the Wigner 
function factorizes and when the state \ip% n } = |l+,0— ) 
is injected, it reads: 



W\ 1+ , -)(a,p,t) = w + (&)W-(i3) 



1 - 4 |a| 2 j exp [-2 |a| 2 l exp [-2 |/3| 2 



(44) 



where a(t) — acosh(<7) — a* sinh(g) and (3(t) = 
/3cosh(g) — f3* smh(g). The variables (^a, (3^j are re- 
lated to (a, 0) through a simple rotation into the quadri- 
dimensional phase space: 

cosh (g) a — sinh (g) (3* — cosh (g) a — sinh (<?) 5(45) 
cosh (g) (3 — sinh (g) a* = cosh (<?) (3 — sinh (<?) f3\46) 

We note that these rotations take the form of the Bogo- 
lioubov transformations which express the time evolution 
of the field operators in the collinear Optical Parametric 
Amplifier. 
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SECTIONS OF WIGNER FUNCTION FOR TWO MODES QI-OPA 
INJECTED QUBIT AMPLIFIED FIELD 
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FIG. 5: Sections of the Wigner function for the injection of 
the single photon state and for the corresponding amplified 
state obtained with a non-collinear amplifier with a gain value 
g = 4. 



C. Measurement of the quadratures with a double 
homodyne 

Quadratures measurement can be obtained by homo- 
dyne technique [25[ largely adopted in the context of 
quantum optics. Quantum fields showing; quad ratures 
entanglement have been realized [H, H3, l38l. l39l liol lllj 
and adopted in several quantum information protocols. 
We now discuss the realization of a double homodyne ex- 
periment in order to investigate the Wigner distributions 
in the QIOPA case. We define the quadrature operators 
in a general form, introducing the phase dependence from 
the variable 8: 



-18 



(47) 



Let us briefly review the scheme of an homodyne exper- 
iment. The impinging field under investigation E\(r,t) 
is combined into a beam splitter with a local oscillator 
EL.o{ r ,t)i usually prepared in a coherent state with 
(3 = |/3|e* e . When f3 is larger than the amplitude of the 
field Ei , the difference between the photon-counts of the 
two output modes of the beam-splitter is proportional to: 



(c^) - (aid) « y/2T{l-T)(Y e ) 



(48) 



where r is the BS transmittivity. By varying the local os- 
cillator phase by tt/2 it is possible to select the measured 
quadrature. 

In the context of the QI-OPA we need to generalize the 
homodyne measurement for the two polarization modes. 
A BS is inserted in the experimental scheme, following 
two different solutions. In the first case (FigJ5]-(a)) at the 



exit of the QI-OPA a couple of waveplates A/4 + A/2 and 
a PBS divide the two orthogonal polarizations which are 
combined with two equally polarized local oscillators. 



(a) 



Local Oscillator 



l«hn> 



QI-OPA 



W4 PBS 



Local Oscillator 



(b) 



l<t>,N> 



QI-OPA 




BS 



IZHHHZr-D,. 



W4 PBS 



Local 
Oscillator 



FIG. 6: Two equivalent schemes for the double homodyne 
measurements. In both schemes, the choice of the basis for the 
analysis depends from the waveplates A/4 + A/2 settings, (a) 
A PBS divides the two polarization components, each compo- 
nent is combined with a coherent radiation (with same polar- 
ization) on a BS. At the output of the PBS the field is revealed 
by two photodiodes. The signals are analyzed following the 
theory of the homodyne technique, (b) A BS combines the 
field with a local oscillator. The two output field impinge on 
a PBS which divides the two polarization components. The 
choice of the analysis basis depends on the waveplates setting 
A/4 + A/2. Signals relative to the different polarization com- 
ponents are detected through two photodiodes and the results 
are pair-correlated according to homodyne theory. 

In the second case (Fig|S]-(b)) the field at the exit of 
the QI-OPA is combined with a local oscillator onto a 
BS; the two output modes are analyzed in polarization 
through a couple of waveplates A/4 + A/2 (same setting 
on both modes) and a PBS and finally detected by a pair 
of photodiodes. The local oscillator polarization must be 
intermediate with respect to the analysis basis, for exam- 
ple if the analysis basis is the linear ±45° rotated one, 
the polarization of the local oscillator can be horizontal 
or vertical. The intensity measurements on both arms 
must be correlated, coupling equal polarizations in or- 
der to obtain the result of (|48|) . The results obtained in 
the two configurations are equal, and the two schemes are 
completely equivalent for the characterization of QI-OPA 
amplified states. 
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D. Optical parametric amplification of N> 1 Fock 
states 

Let us now consider the two-photon input state \ipi n ) = 
|2+,0— ) the characteristic function reads: 

X aAv,t,t) = ex P [-i(h(*)| 2 +K(*)| 2 )] fl - W*)l 2 + lj ^r~ 



and the Wigner function becomes: 



(49) 



|o) 

BS 



Detector 



R+T=l 



FIG. 7: Schematization of the decoherence model by a linear 
beam-splitter of transmittivity T. 



W\ 2 +fl-)(a,/3,t) 



1-2IA, 



-|A| 2 



l\A A + A B f 



(50) 



A three dimensional plot of this function is obtained (sec 
Fig. 4). 

When both polarization modes are one photon injected 
\i>in) — 1 1 — I— , 1 — ) , the characteristic function reads: 



xmOj.£.*) = 1 



\v(t)\ 



and Wigner function is: 



l€(*)l 



(51) 



W {1+<1 _ } (a,p,t) = - [l-2|A| 2 4-|A A + A B | 2 x 



x (1-|A A 4-A s | 2 )le- |A| 



(52) 



As in the single mode case, we can generalize the re- 
sults obtained for a generic Fock-state as input \ipin} — 
\N+,M— ). The Wigner functions of the amplified field 
with these generalized seeds read: 



(-1) 



M+N 



L 



N 



Ar| 2 ) 



L M (|A B -A A | 2 ) e- 2 ^ 2 



. (53) 

The Wigner function approach given in this Section to 
the QI-OPA device allows to stress the quantum proper- 
ties of the amplified field. The negativity of the Wigner 
function can be deduced by the explicit general ex- 
pression of Eq. (|53p . Indeed, in analogy with the un- 
amplificd Fock-states \N+,M—) [22j], the Wigner func- 
tion W\ n+. m -\ (a, 0, t) is the product of two Laguerre 
polynomials [42j Lm and Lm, which are negative in sev- 
eral region. Furthermore, the amplified states possess a 
high degree of squeezing in the field {a, 0\ quadratures, 
at variance with the injected states. 



photon states are smoothed and cancelled when a deco- 
herence process, i.e. a "system-enviroment" interaction, 
is affecting their time evolution. More specifically, in the 
specific case of optical fields, the main decoherence pro- 
cess can be identified with the presence of lossy elements, 
as for example photo-detectors. Such process is mathe- 
matically described by an artificial beam-splitter (BS) 
scattering model, since this optical element expresses the 
coupling between the transmission channel and a differ- 
ent spatial mode. The same analysis is carried out also 
on a different class of MQS based on coherent \a) states, 
to emphasize analogies and differences. 



A. The lossy channel model 

As said, losses are analyzed through the effect of a 
generic linear beam-splitter (BS) with transmittivity T 
and reflectivity R = 1 — T, acting on a generic quantum 
state associated with a single mode beam: FigfT] [H, H3, 
EE Ell- The action of the lossy channel on a generic 
density matrix p is obtained by the application of the 
BS unitary transformation and by the evaluation of the 
partial trace on the BS reflected mode (R-trace). The 
linear map describing the interaction is expressed by the 
following expansion [44], [45| : 



where the M p operators are: 



M p = Rvl^'^l 2 - 



(54) 



(55) 



To evaluate the average values of operators, the lossy 
channel can be expressed in the Heisenberg picture, ex- 
ploiting the BS unitary transformations and performing 
the average on the initial state p. 



IV. FOCK-SPACE ANALYSIS OF 
MACROSCOPIC QUANTUM SUPERPOSITIONS B. De-coherence on the quantum superposition of 
OVER A LOSSY CHANNEL coherent \a) states 



In the present Section we analyze how the peculiar 
quantum interference properties of the amplified single- 



In this Section we investigate the evolution over a lossy 
channel of the quantum superpositions of coherent states 
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(CSS) [ag: 



|^>=^r-±=(K 



± \ae~ l f)) 



(56) 



with 



real 



and 



J v 



/]_ ±e -2|a| sin y cos [| a |2 sin 2 ¥ ,]J " is an appro- 
priate normalization factor. The two states with 
opposite relative phases |$+) and |\&^) are orthogonal 
when |a| 2 sin 2 t^ > 1 (FiglSJ) . In such case the two 
components \ae lv ) and \ae~ lv ) are distinguishable. This 
class of Macroscopic Quantum Superposition presents 
several peculiar properties, such as squeezing and 
sub-poissonian statistics, which can not be explained by 
the characteristics of the coherent la) states. 



sin tp 




FIG. 8: Schematic representation in the bidimensional phase- 



space {X, Y} of the coherent states 



') which are the 



components of the MQS treated in this section. The cir- 
cles represent the two-dimensional projection of the Gaus- 
sian functions which represent the Wigner function of these 
states. The distance between the center of the two Gaussian 



is cL 



: 4|q 



2 ■ 2 

sm tp. 



We now proceed with the analysis of the density ma- 
trix after the propagation over the lossy channel. In 
the following we assume \a\ 2 sin 2 tp > 1, hence J\f^ ~ 1. 
The density matrix of the quantum state after the beam- 
splitter transformation is: 



R=0.1 




FIG. 9: (a)-(d): Plot of the distribution of the number of 

photons in the \^?) state for a — 4, corresponding to an 
average number of photons (n) = 16, for reflectivities R = 
(fig.2-a), R = 0.1 (fig.2-b), R = 0.5 (fig.2-c) and R = 0.8 
(fig-2-d) 



the density matrix after losses reads: 

Pc = \ </^| + \Pe~ llp ) (Pe- iV \ + 

_j_ e -2_R|a| 2 sin 2 Ve i_R|a| 2 sin2^|^ e J^^ /^g-tp|_|_ (58) 
_|_ e -2i?.|a| 2 sin 2 tp £ -iR\a\ 2 sin 2tp |^g— ^9e tp |J 

Let us analyze first the case tp = The distribution in 
the Fock space exhibits only elements with an even num- 
ber of photons for |^+) or an odd number of photons for 
| , 5_). This peculiar comb structure is indeed very frag- 
ile under the effect of losses since the R-trace operation 
must be carried out in the space of the non- orthogonal 
coherent-states. In fig |9] are reported the photon number 
distribution for different values of R with an initial av- 
erage number of photon equal to (n) = 16 and tp = \. 
We observe that for a reflectivity R = 0.1, corresponding 
to about ~ 1.5 photon lost in average, the distribution is 
similar to the Poisson distribution associated to the co- 
herent states. Such graphical analysis points out how the 
ortoghonality between \^f+) and ) quickly decrease as 
soon as R differs from 0, since the phase relation between 
the components \a) and | — a) becomes undefined. 



ff ±V = \ (\Pe^)cc{Pe^\ ® \i 1 e^) DD {t 1 e^\ + 

+ \Pe-^) CC {0e-^\ ® \i-ye-^) D D {ije~^\ + 
± W lfi ) CC (f3e-^\ <g> \ije^) D D (ije'^\ + 
± IPe-^c c (I3e iv | ® \tje~^) D D (vye^\) 

(57) 

with (3 — y/Ta and 7 = \/i?a. The output state over 
the transmitted field is obtained by tracing the density 

matrix p 1 over the mode d. The final expression for 



C. De-coherence in a lossy channel for equatorial 
amplified qubits 

Analogously to the previous Section, in which the effect 
of losses where analyzed for the MQS of coherent state, 
we begin our treatment of QIOPA amplified states with 
the evaluation of the density matrix after the propagation 
over a lossy channel. 

Before the lossy process, the density matrix of the state 
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pv = has the form: 

OC 

P V = E 7y7fc g l(2* + lV, (2j> ± ) ((2fc+%, (2q)tp x \ 

i,j,k,q—0 

. (59) 

We note from this expression that only elements with an 
odd number of photons in the n v and an even number in 
the if v± polarization are present. Furthermore, in Fig llOt - 
(a) we note that the photon number distribution presents 
a strong unbalancement due to the quantum injection 



of the 7?^ single photon. Indeed, the QIOPA seeded by 
a photon with equatorial polarization acts as a phase- 
covariant optimal cloning machine, and is stimulated to 
generate an output field containing more photons in the 
polarization of the injected seed. Let us now analyze 
the effects of the transmission in a lossy channel for the 
equatorial amplified qubits by plotting the photon num- 
ber distributions. The output density matrix after the 
transmission over the lossy channel is the sum of four 
terms: 



Pt 



= E g^W^ 1 )^ 2 ^ 

i,j,k,q=0 (60) 

+ (%) ijkq |(2* + IV, (2j + l)ip ± ) ((2k + IV, (2q + + (p*) ijkq |(2*>, (2j + l)<p ± ) ((2k)<p, (2q + l)<p ± \ 



r 



The details on the calculation and on the expressions of 
the coefficients are reported in appendix [A"l 

Let us now analyze this result. When the original 
state propagates through a lossy channel, the first ef- 
fect at low values of R is the cancellation of the peculiar 
comb structure (FigllOIfa)) given by the presence in the 
density matrix (|59p only of terms with a specific parity 
|(2* + l)(f, (2j)(p±) ((2k + l)cp, (2q)<p±\, similarly to the 
coherent state MQS's previously studied. However, at 
progressively higher values of R, the distributions in the 
Fock space remain unbalanced in the polarization of the 
injected photon (FigfTOl-(a)). The resilience of this un- 
balancement allows to distinguish the orthogonal macro- 
qubits {l^), l^^^)} even after the propagation over the 
lossy channel, by exploiting this property with a suitable 
detection scheme, such as the OF device reported in [2l[ . 
All these considerations will be discussed and quantified 
later in the paper in Section IVIII by analyzing the dis- 
tinguishability of such states as a function of the lossy 
channel efficiency T. 



D. De-coherence in a lossy channel for amplified 
tth,t(v qubits 

For the sake of completeness, in this Section we shall 
analyze the evolution of |«& ) and \& v ) amplified states. 
As a first remark, we note that the collinear Optical Para- 
metric Amplifier is not an optimal doner for states with 
tth and iry polarization, and the output states do not 
possess the same peculiar properties obtained with an 
equatorial injected qubit. The density matrix of the ) 
amplified state is: 



\<$> H ) (§ H \ 



1 



E r ,l+m V^TT\/^TT 



\(n + \)H,riV) ((m + l)H,mV\ 



In FigfTUKb) we plotted the photon number distribution 
of this state (R=0). We note that the tth amplified state 
does not possess the same unbalancement of the equato- 
rial macro-qubits analyzed in the previous section. 



After the propagation over the lossy channel, the den- 
sity matrix reads: 



Pt 



EEEC 

i=l j=0 k=0 



p?: 



ijk 



\iH,jV) (kH,(k + j-i)V\+ 



+EEE 

i—0 j—i k—0 



Pt] 



tjk 



\iH,jV) (kH,(k + j-i)V\ 



(62) 

where details on the calculation and on the coefficients 
are reported in appendix [Bj The effect of the propaga- 
tion over the lossy channel is shown in Fig |10l -(bV The 
original distribution for R = is pseudo-diagonal, corre- 
sponding to the presence only of terms \(n + l)H,nV). 
Here the difference of one photon between the two po- 
larization is due to the injection of the seed. For val- 
ues of R different from 0, the distribution is no longer 
pseudo-diagonal and this characteristic becomes progres- 
sively smoothed. Furthermore, the absence of the un- 
balancement in the photon number distribution typical 
of the equatorial macro-qubits does not allow to exploit 
this feature to discriminate among the orthogonal states 
We then expect that these couple of states 
possess a lower resilience to losses than the equatorial 
macro-states analyzed in the previous section. 



(61) 
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FIG. 10: (a) Probability distribution in the Fock space (n^,n^) for the amplified |$^} state of a generic equatorial qubit for 
different values of the transmittivity. (b) Probability distribution in the Fock space (nn,nv) for the amplified \& H ) state for 
different values of the transmittivity. All distributions refer to a gain value of g = 1.5, corresponding to an average number of 
photons (n) w 19. 



V. WIGNER FUNCTION REPRESENTATION 
OF COHERENT STATES MQS IN PRESENCE OF 
DE-COHERENCE 



For the sake of clarity, we briefly review previous re- 
sults on the Wigner functions associated to coherent su- 
perposition states after the propagation over a lossy chan- 
nel. We start from the general definition of the Wigner 
function for mixed states 1231] : 



Considering the density matrix of the CSS after losses 
we obtain: 



(W We , v) (X, Y) + W\ Pa -^){X, Y)+ 

rint 



±W™ ± (X,Y) 
Pc 



(64) 

In the last expression, the first two components are anal- 
ogous to the diagonal ones of the unperturbed case [IH , 
and can be written as : 

W We± , v) (X, Y) = I e -(*-v^)V( y TvW,) 2 (65) 



i r°° 

Wp{X, Y) = - d^{X - £\p\X + (63) 

J-oo 



where X 2 ^ — 2|a| 2 cos 2 </j and Y — 2|a| 2 sin 2 ip. Hence 
losses reduce the average value of the quadratures X and 
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Y. 



The interference contribution reads 
2 

(A. ) i 



— RY 



X COS 



2V2aVT sin w ( X - Q ( 27 ^ 1) cos „, ) 
V V2T V 



which is strongly reduced in amplitude by a factor pro- 

portional to e . 

In FigQj] are plotted the Wigner functions associated 
to different values of i?, for the same initial conditions 

= § an d a = 6. As expected, by increasing the degree 




FIG. 11: Plots of the Wigner function for a — 6 and 

for different values of the riflectivity (a) R = 0, (b) R = 0.01, 

(c) R = 0.02 and (d) R = 0.1. 

of losses the central peak is progressively attenuated up to 
a complete deletion of the quantum features associated to 
the negativity of the Wigner functions. We observe that 
the damping factor e ~ 2R \ a \ Sln v of the coherence terms 
derives from the exponential decrease of the non-diagonal 
terms of the density matrix (|58[) . 

More specifically, we now focus on the = § case, 
i.e. the \a) ± | — a) state. In Fig. [12] we report the 
plots of the Y = section of the Wigner function for 
different values of the reflectivity R, which corresponds 
to the interference pattern. 

At low R, the amplitude of the oscillation is exponen- 
tially damped. This exponential factor e~ 2fl ' Q ' sm v is 
responsible for the fast decrease in the negativity of the 
Wigner function, but the alternance of positive and nega- 
tive regions is maintained. However, when R approaches 
the 0.5 value, the interference pattern is progressively 
shifted towards positive values in all the X-axis range, 
and at R = 0.5 it ceases to be non-positive. This evolu- 
tion depends on the presence of the W\p e ±%v\{X, Y) di- 
agonal terms, which are not exponentially damped in 
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FIG. 12: Y = sections of the Wigner function for the quan- 
tum superposition of coherent states, for a = 6 and tp = 2 
and different values of R. (a) R = 0: Unperturbed interfer- 
ence pattern, (b) R = 0.05: Exponentially damped interfer- 
ence pattern, (c) R = 0.5: Transition to a completely positive 
Wigner function, as the interference pattern is shifted towards 
positive values, (d) R = 0.7: The interference pattern in the 
central region is deleted by decoherence. 



amplitude and for R ~ 0.5 become comparable to the 
interference term. This transition is graphically shown 
in Fig |131 where the negativity is plotted as a function 
of R. This quantity has been evaluated by calculat- 
ing the value of the Wigner function in the first mini- 
mum of the cosine term of W!" ± (X, Y) . We note the 

Pa _ 

transition from negative to positive at R = 0.5, where 
for higher values this point ceases to be the minimum 
of the complete Wigner function in the {X, Y} plane. 
This property can be derived explicitly by the complete 
form of the Wigner function. For tp — ^, Y" = and 
Xq = - ^ ^ — g ^- , corresponding to the minimum of the 



cosine cos 



2V2a 



Tsintp [X - 



q(2T-l) 
V2T 



cos tp 



MI 



W f + (X ,0) = ^e 



-2(1-«)H' 



-2|q| 2 (1-H) 



we get: 



-2\a\ 2 R 



< if R < 



> if R > f 

(67) 

thus obtaining the desired result. 

As a general statement, it is well known that the neg- 
ativity of the " quasi-probability" Wigner function of a 
state p is a sufficient, albeit non necessary, condition for 
the " quantumness" of p 22]. 



VI. WIGNER FUNCTION REPRESENTATION 
OF THE PHASE-COVARIANT QUANTUM 
CLONING PROCESS IN PRESENCE OF 
DE-COHERENCE 

Combining the approach of Section III and the lossy 
channel method introduced in Section IV, we derive the 
analytical expressions for the Wigner function of the QI- 
OPA amplified states in presence of losses. The calcula- 
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FIG. 13: Transition from non-positive to completely positive 
Wigner function as a function of the reflectivity R, for a num- 
ber of average photons of (n) ~ 6.5 (a) and (n) « 53.5 (b). 
The negativity has been evaluated as described in the text. 
We note that, while the interval of R in which the Wigner 
function is non positive is independent of the number of pho- 
tons, the amount of negativity decreases exponentially with 
(n) due to the exponential factor e~ 2fl ' Q ' sm v in the inter- 
ference term W%*±{X,Y). 
p c 



tion will be performed in the Heisenberg picture, starting 
from the evaluation of the characteristic function of the 
BS-transmitted field. This analysis will be performed for 
both the single-mode degenerate amplifier, i.e. a single- 
mode squeezing Hamiltonian, and for the two-mode op- 
tical parametric amplifier, in which the polarization de- 
gree of freedom plays an important role as stressed in the 
Fock-space analysis previuosly performed. Finally, we 
shall focus our attention on the negativity of the Wigner 
function, evidence of the non-classical properties of this 
class of states. 



A. Single-mode degenerate amplifier 

Let us first analyze the case of a single mode- 
degenerate amplifier. We begin our analysis by evalu- 
ating the characteristic function in presence of a lossy 
channel. The operators describing the output field can 
be written in the Heisenberg picture in the form: 



(68) 



where aJ(t) = d^cosho + ashing is the time evolution 
of the field operator in the amplifier, and R = 1 — T. 
Hence, the characteristic function for a generic input 
state l-ZV^lC^s can be calculated as: 

XN (v,R,t) = (^0| e " £t M-"* £ (*>|iV,0) (69) 
Inserting the explicit expression |68|) for the output field 
operators, we obtain the expression: 



X N(r ] ,R,t)=e-^W 2 e-i R M\ 



(70) 



In this expression, the transformation (77,77*) — > 
(r](t) , T]* (t)) has the form r)(t) = r\C — rj*S, equivalent 
to the ideal case previously analyzed. The Wigner func- 
tion of the output field is then obtained from its defi- 
nition (|20[) as the two-dimensional Fourier transform of 
the characteristic function. Inserting the explicit expres- 
sion of xn(ti, R,t), and by changing the variables with 
the transformation r/ — > 77(f) with unitary Jacobian we 
obtain: 



W\ N )(a,R,t) 



1 N 

^2 



n=0 



d 2 v(t)\v(t)\ 



2n-e\n(t)\ 2 



xexp[- K (77 2 (f)+7f 2 (<))] 
x exp [— a* r)(t) + arf(t)} 



where: 



N\ (-1)™ 



rj-vn 



1 + 2RS 2 



k = -RCS 

a = aC — a 



"S 



(71) 

(72) 

(73) 

(74) 
(75) 



In order to evaluate the integral in the previous expres- 
sion, we recall the following identity [33L l47l| : 



7T 



y/T 2 - ±\1V 



exp 



-fia — va — z a + za 

v 2 1 ,,_*2 1 ^1,12" 



fiz 2 + vz* 2 + t\z\ a 

T 2 — AflV 



(76) 



From the latter equation, we can derive the following 
useful identities: 
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2 * 2 * i * 

-\ia — va — z a + za 



I n (m, v,t,z) = ^- J d 2 a |a| 2 " e - T|Q|2 exp [- 

= -(-1)" 8 " I d 2 ae~ T ^ 2 exp \-p,a 2 -va* 2 - z*a + za* 
it v ' dz n dz* n J F L p 

/xz 2 + i/z* 2 + t|z| 2 



= (-1)" 



<9 2ri I 1 

r 



4/w 



(77) 



These expressions can then be used to evaluate the 
Fourier transform of the characteristic function (|70|) . 
leading to the Wigner function: 



1 N 

W lN) (a,R,t) = -Y / ~ h 

n=0 



c n I n (k, k, e, a) 



(78) 



As a first case, we analyze the evolution of the Wigner 
function of a squeezed vacuum state, corresponding to 
the value N = in the previous expression. Analogously 
to the unperturbed case, the quadrature variables for the 
single-mode OPA are defined by a = X+zY. The Wigner 
function for the squeezed vacuum then reads: 



W\ 0) (A, Y, R,t) = - 1 

1 ' T ^1 + 4(1 -R)RS 2 

(X 2 



x exp 



-2- 



< 2,1 + Y 2 e 29 ) + 2RS (X 2 e 
1 + 4(1 - R)RS 2 



Y 2 e s) 



(79) 



From this quantity, we can explicitly calculate how the 
degree of squeezing changes with an increasing value of 
the parameter R. We obtain for the fluctuations of the 

two field quadrature operators |a, Y~j (Fig ITU- (a)): 




FIG. 14: (a)-(b) Uncertainty relations for the quadrature op- 
erators in the spontaneous emission case for a single mode 
OPA. Left figure: Fluctuation on the X quadrature (red 
straight line) and on the Y quadrature (green dashedline) as 
a function of the parameter R of the lossy channel. Right 
figure: Uncertainty relation AX AY as a function of the pa- 
rameter R. (c)-(d) Uncertainty relations for the quadrature 
operators in the single-photon amplified case for a single mode 
OPA. Left figure: Fluctuation on the X quadrature (red 
straight line) and on the Y quadrature (green dashed line) as 
a function of the parameter R of the lossy channel. Right 
figure: Uncertainty relation AX AY as a function of the pa- 
rameter R. All figures refer to a gain value of g — 3. 



AX = -y/Te 2 a + (1 - T) 



AY 



-Vfe=*9 + (1-T) 



(80) 
(81) 



This two operators do not satisfy anymore the minimum 
uncertainty relation, as in the unperturbed case, after 
the propagation over the lossy channel: Fig[T4lfb). This 
is due to the additional poissonian noise belonging to the 
photon loss process. 

A similar behaviour is obtained for the single-photon 
squeezed state. The Wigner function for N = 1 reads: 



W {1) (X,Y,R,t 



1 



P m (X,Y,R,t) 



x exp 



7T y/1 + 4(1 - R)RS 2 
(X 2 e- 2a + Y 2 e 2 3) + 2RS (X 2 eS + Y 2 ^) 



1 + 4(1 - R)RS 2 



(82) 



where: 

P w (X,Y,R,t) = l 



4(1 -R) 



(X 2 e~ 2 ° 
(X 2 e- 2 9 



Y 2 e 29 ) 
- Y 2 e 2 s ) 



1 + 4(1 - R)RS 2 
+ 2(1 + 2RS 2 ) x 



'l + 2i?S' 2 



2RS (A 2 



Y V) 



1 + 4(1 - R)RS 2 



(83) 



The region where the Wigner function is negative be- 
comes smaller when the parameter R of the lossy chan- 
nel is increased. In Fig[TS] we report the plots of 
W\ i) (X, Y, R, t) for different values of the reflectivity R. 
As a first effect, the negative region is deleted for a re- 
flectivity R = 1/2: FigElKd). Then, the form of the 
distribution remain unchanged until the reflectivity be- 
comes close to 1 and all the photons present in the states 
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FIG. 15: Wigner function of a single-photon amplified state in 
a single-mode degenerate OPA for g — 3. (a) (R=0) Unper- 
turbed case, (b) (R=0.005) For small reflectivity, the Wigner 
function remains negative in the central region, (c) (R=0.1) 
The Wigner function progressively evolve in a positive func- 
tion in all the phase-space, (d) (R=0.5) Transition from a 
non-positive to a completely positive Wigner function. 



are lost: R(n) ~ (n). As a further analysis, let us con- 
sider the value at X = and Y — 0, in which the Wigner 
function has the maximum negativity. We obtain that 
W w (0,0, R,t) < for R < §, showing that the nega- 
tivity is maintained in that range of the lossy channel 
efficiency. This property can be analyzed by the two- 
dimensional contour plots of figfTHl 

As a further analysis, starting from the Wigner func- 
tion, we can evaluate the fluctuations in the quadrature 
operators after losses (FigQ3r(c): 



AX = \pJZTe 2 a + (1 - T) 
AY = i v /3Te- 2 f + (l-T) 



(84) 
(85) 



As in the squeezed vacuum case, poissonian noise added 
by the lossy channel increases the fluctuations AX AY 
for any non-zero value of the losses parameter R: FigQjJ- 
(d). ' 
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FIG. 16: Contour plots of the Wigner function of a single- 
photon amplified state in a single-mode degenerate OPA for 
g = 3. (a) (R=0) Unperturbed case, (b) (R=0.005) For 
small reflectivity, the Wigner distribution begin to evolve in a 
double-peaked function with a negative central region, anal- 
ogously to the CSS case, (c) (R=0.5) Transition from the 
regime with a negative central region and a completely posi- 
tive distribution (d) (R=l) Evolution of the Wigner function 
into a vacuum state when all photons are lost. 



cess and to the propagation in the lossy channel, takes 
the form: 



ff H (t) = Vl - Ra ] H (i)+i 
c ] v (t) = Vl - Ra\,(t)+i 



Rb 



H 



(87) 



Rb. 



where the expressions of the equation of motion for 
(t) are ([4]). We now proceed by writing the char- 



l H,V 



acteristic function in the {7? + ,7?_} basis and by substi- 
tuting the relations (|87I88[) for the time evolution in the 
Heisenberg picture of the field operators. The character- 
istic function then reads: 



B. Two- mode collinear amplifier 

The previous calculation, performed in the case of a 
single-mode degenerate OPA, can be again used to ana- 
lyze the collinear QI-OPA case. The characteristic func- 
tion with a generic input two-mode state \^i n ) can be 
evaluated as: 



XN, M (ri,i,R,t) = (* m |e"4W->*«(*) ( 



(86) 



Analogously to the single-mode OPA case, the time evo- 
lution of the field operators, due to the amplification pro- 



XN,MKn,£,,R,t) = (*m|e 

ivO^-K) j;t -iVRjy'+H"'-) ft 



V2 + V2 + 



V2 



73 "+ 



e V2 73 



V2 



(89) 



The transformation between (rj, rf , £, £*) — > 
(ri(t),rf(t), £,(£),(,* (t)) is the same (life]) of thc 
unperturbed case. The evaluation of the average on an 
input injected state \N+, M— ) in the QI-OPA and on 
the vacuum-injected port b of the beam-splitter leads to 
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the following result: tion. Analogously to the single-mode case, we proceed 

with the calculation by evaluating the Fourier integral 
Xn M(f?>£j-ft>*) = e _ ' it (' , ' + ' 5 ' 2+ ^~''' 2 ) that, after changing the integration variables with the 

N M 

transformation (ij^.^f) -» (*?(*),»/*(*),€(*),€*(*)), 
„Af can be written as: 



{N M 
E E c « c " 
n=0 m=0 



(90) 



The Wigner function is then calculated as the 4- 
dimensional Fourier transform of the characteristic func- 



W\ N+! M-}(a,/3,R,t) 



17^ 
N 



d 2 V (t) I d^(t)|J|exp 

M 



1 



E c * 

Vn=0 



2n 



E 



2 

K(*)i 2m 



1 



exp 



:«(^(t)+„*2(i)-^(*)-r a (*)) 



exp 



V2 



(T7*(t)a-tj(t)5*) 



exp 



^ (f(^-?(t)/?' 



(91) 



where the parameters e, k have been defined in (|73H74[) . 
and the transformation between a — > a and /3 — > /3 is: 



a 



ft 



1 

71 
i 

71 



[{a + ft)C-(a* +ft*)S] (92) 
[G8-a)C+(/T-a*)S] (93) 



The derived integral relations ([77|) lead to the final result: 

W\ff +t M-) I 



! / w _ \ 

, (a, /3, R, t) = -2 I ^ CrJniK, K, £, a) J 
n Vn=0 / 



/ M \ 
X I ^ ] C m I m ( — K,—K,£,f3) 
\m=0 



(94) 

Let us analyze the spontaneous emission case, when 
N = M = 0. Adopting the same definition of the phase 
space used in the unperturbed case in Section IIII Bl we 
obtain: 



-R) 



x exp 



x exp 



(l + 2RS 2 + 2RCS) X 
\ + AS 2 R(l- R) 

(1 + 2i?S 2 - 2RCS) Y 2 e 2 9 
1 + AS 2 R(1-R) 



(95) 



As in the single mode case, the degree of squeezing in 
the quadrature operators is progressively decreased by 
the poissonian noise introduced by the lossy channel. 



Furthermore, let us analyze the case of the single pho- 
ton amplified states, i.e. N = 1, M = and N = 0, M = 
1. The Wigner functions after losses for this quantum 
states are: 



W\i+, -){X,Y,R,t) 
W|o+,i-) (A, Y, R, t) 



WV,o-)^|i+,o-}(A,i?,ii96) 
W|o + ,o->-fV,i->(^*)(97) 



where P| 1+i0 -) (A, R, t) and P\ 0+ ^_)(Y,R,t) are the fol- 
lowing 2 nd order polynomials: 



Pu +fi _){X,R,t) = 1 



4(1 - R) 



1 



1 + 4S' 2 i?(l - R) 
-X 2 e- 2g ) - RS 2 + (1 
2RS 2 + 2RCSX 2 e~ 29 ' 



2RS 2 



1 + 4R(1 - R)S 2 



(98) 



4(1-*) 



1 + 45' 2 i?(l - i?) 



(l + FV 9 ) - RS 2 + (l + 2i?S" 2 



1 + 2i?S 2 - 2RCSY 2 e 2 9 



1 + 4R(1 - i?)^ 2 



(99) 



In FigQ~7] we report the Wigner function 
W\i+o-\(X,Y,R,t) for different values of the re- 
flectivity R. The evolution of this distribution is similar 
to the single-mode OPA case analyzed previously. 
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FIG. 17: Wigner function of a single-photon |1+) amplified 
state in a two-mode degenerate OPA for g — 3. (a) (R=0) 
Unperturbed case, (b) (R=0.005) For small reflectivity, the 
Wigner function remains negative in the central region, (c) 
(R=0.1) The Wigner function progressively evolve in a posi- 
tive function in all the phase-space, (d) (R=0.5) Transition 
from a non-positive to a completely positive Wigner function. 



Resilience of quantum properties after 
decoherence 



The Wigner functions calculated in the previous sec- 
tion allows to obtain a complete overview of the phase- 
space properties of the QI-OPA amplified states after the 
propagation over a lossy channel. In particular, we focus 
on the negativity of the W-representation in the specific 
case of a single-photon injected qubit. As said, the pres- 
ence of a negative region in the phase-space domain is one 
possible parameter to recognize the non-classical proper- 
ties of a generic quantum state. 

Let us now consider the expression ([95)) for the polyno- 
mial P|i_|_ ) o-\(X, R, t). In the lossless regime considered 
in Fig[4l the Wigner function takes its minimum value 
in the origin of the phase-space (X=0,Y=0). In presence 
of losses, the Wigner function remains negative in the 
origin for R < i. This behaviour is shown in the plots 
of the Y = sections of Wji +i g_\ (X, Y, R, t) reported in 
FigfTSl where for the value of R = 5 ceases to be nega- 
tive. This is also shown in FigfTSl that reports the value 
of l / F|i +j o-) (0, 0, R, t) as a function of the reflectivity R of 
the beam-splitter that models the lossy channel. Finally, 
this behaviour, as for the quantum superposition of co- 
herent states in the previous Section, can be analytically 
obtained by calculating Pi 1+>0 _\(X, R,t) in X = 0. We 
obtain: 



P\i+,o-)(X, R,t) — 



2R- 1 



1 + 4R(1 - R)S 2 



< 
> 



if 



R< I 
if R> I 
(100) 

This analysis explicitly shows the persistence of quantum 
properties in the QI-OPA single-photon amplified state 
even in presence of losses. However, the negativity of the 
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FIG. 18: Y = section of the Wigner function for QIOPA 
single photon amplified states for (n) fa 36. (a) R = 0: Unper- 
turbed Wigner function, negative in the origin, (b) R = 0.05: 
The negativity in the origin progressively decreases due to the 
coupling with the enviroment. (c) R — 0.5: Transition from 
a non-positive to a completely positive Wigner function, (d) 
R = 0.7: The Wigner function is positive in all the {X, Y} 
plane. 
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FIG. 19: Transition from a non-positive distribution to a com- 
pletely positive Wigner function for QI-OPA single-photon 
amplified states in presence of decoherence, as a function of 
the beam-splitter reflectivity R. Average photon number of 
(a) (n) w 6.5 and (b) (n) w 53.5. The negativity is evaluated 
as the value of the Wigner function in the origin. 



Wigner function is not the only parameter that reveals 
the quantum behaviour in a given system. Hence, more 
detailed analysis have to be performed including different 
approaches in order to investigate the R > | regime. 

We conclude our analysis by comparing the behaviour 
of the two different macroscopic quantum superpositions 
analyzed in this paper. In Fig. [201 we compare the de- 
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FIG. 20: Absolute value of the negativity, evaluated as the 
minimum of the Wigner function, in the reflectivity range 
< 7? < |. The blue straight lines correspond to QI-OPA 
amplified states, while red dashed lines to the \a) ± | — a) 
quantum superposition, (a) (n) ~ 6.5 and (b) (n) ~ 53.5. 



crease of the Wigner function negativity for the QI-OPA 
single-photon amplified states and for the superposition 
of coherent states, for two different values of the aver- 
age number of photons. We observe that the QI-OPA 
solution possesses a higher resilience to losses, i.e. a 
slower decrease of the negative part, with respect to the 
\a) ± | — a) states. However, the Wigner function for both 
quantum superposition ceases to be negative at R = 4. 



VII. PERSISTENCE OF COHERENCE IN 
MACROSCOPIC QUANTUM SUPERPOSITIONS: 
FOCK-SPACE ANALYSIS 

In this Section we perform a complementary analysis 
in the Fock-space of the macroscopic quantum superpo- 
sition generated by the quantum cloning of single-photon 
states, by applying two criteria [29| based on the concept 
of distance in Hilbert spaces, more specifically the Burcs 
metric. This approach allows to quantify from a differ- 
ent point of view the different resilience to losses that this 
quantum superposition posses in contrast to the fragility 
of the | a) states MQS. 

First, we introduce the coherence criteria and discuss 
their interpretation. Then we apply this last approach 
to the two different MQS under investigation showing 
analogies and differences. Furthermore, with an oppor- 
tune POVM technique, based on the O-Filter device in- 
troduced in [2(| HH, the properties in the Fock-space 
of QIOPA amplified states can be exploited to obtain 



a higher discrimination in the measurement stage, at the 
cost of a lower events rate. 



Criteria for Macroscopic Quantum 
Superpositions 



Metrics in Hilbert spaces. In order to distinguish 
between two different quantum states, we need to define 
a metric distance in the Hilbert space. An useful pa- 
rameter to characterize quantitatively the overlap of two 
quantum states is the fidelity between two generic density 



matrices p and a, defined as: J-(p, <?) = Tr 2 



piap* 



[4c| . This parameter reduces to \<p)) — \(tlj\ip)\ 2 

for pure states, and is an extension of the scalar prod- 
uct between quantum states to the density matrix for- 
malism. We have < T < 1, where J- = 1 for iden- 
tical states, and T — for orthogonal states. This 
quantity is not a metric, but can be adopted to define 
two different useful metrics, which are the angle distance 
Da(p,^) = arccos !F(p, d) [H, and the Bures distance 



D(p,a) = Jl-y/T(p,a) 



(101) 



Furthermore, the fidelity can also be used to calculate a 
lower and an upper bound for the trace distance, defined 
as Dj-(p,5) = iTr \p— a\ and related to the fidelity by 

[l|: l-y/F(p,a) < D(p,a) < ^l-T{p,d). In this Sec- 
tion we will adopt the Bures distance as a metric in the 
quantum state space, as it is connected to the probability 
of obtaining an inconclusive result with a suitable Pos- 
itive Operator Valued Measurement (POVM) [H, 
which is y/F (|0), = \{i/)\<l>)\ for pure states . 

Distinguishability, MQS Visibility. Let us charac- 
terize two macroscopic states \4>i) and \4>2) and the cor- 
responding MQS's: |0±) = ^ (|0O ± |0a)) by adopting 
two criteria. I) The distinguishability between \<fii) and 
|(/> 2 ) can be quantified as D(|</>i), \4>2})- II) The " Visi- 
bility", i.e. "degree of orthogonality" of the MQS's 
is expressed again by: D (|</> + ), \ 4>~))- Indeed, the value 
of the MQS visibility depends exclusively on the relative 
phase of the component states: \4>i) and \<j)2)- Assume 
two orthogonal superpositions l^*): D (\<fi + ), \4>~)) = 1. 
In presence of losses the relative phase between |0i) 
and |02) progressively randomizes and the superpositions 
\4> + ) and \<fi~) approach an identical fully mixed state 
leading to: D(\cf) + }, \(f)~)) — 0. The aim of this Section 
is to study the evolution in a lossy channel of the phase 
decoherence acting on two macroscopic states \4>i) and 
|02) and on the corresponding superpositions l^*) and 
the effect on the size of the corresponding D (\(f>i), (^2)) 
and£>(|^+),|r»- 
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B. Bures distance for \a) states MQS 

In order to quantify the loss of coherence in the macro- 
scopic superposition of coherent states, we estimate their 
relative Bures distance. For the basis a states we have: 

D{\ae zv ), \ae- lip )) = y/l - e - 2T M 2 si " 2 * (102) 

The distinguishability between these two states keeps 
close to 1 up to high values of the beamsplitter reflec- 
tivity, when almost all photons are reflected. 

Let us now estimate the distance between the super- 
position states. We consider the following condition 
|(/3e^|/3e-^)| = e -2T|a| 2 sin 2 v _ Q) which corresponds 
to T\a\ 2 sin 2 tp > 1. Hence, except for very small T, 
the coherent states \(3e ±ttp ) after the propagation over a 
lossy channel remain almost orthogonal. In such situ- 
ation we can associate to {\(3e tv ), \f3e~ lv )} the two or- 
thogonal states of a qubit {|0), |1)}. Let us introduce the 
parameters \ = 2|a| 2 i?sin 2 tp and ip = |a| 2 -Rsin2<y9. The 
density matrices after losses can be represented as 2 x 2 
matrices associated to the qubit state: 

Pc =2\±e-*e-* 1 ) (1 ° 3) 

To estimate the fidelity T (pc V , Pc V ) we nee d to calcu- 
late: 

{P^Pc V {P + c V Y = \^~f X !_° e - 2x ) (104) 
Hence we get: 

F(pc V ,fc' P ) = l-e" 2x = 1 - e - 4 H 2 W^ (105) 
From the definition (jlOip . we found: 

D = V 1 - e" 4J? l Q l 2 sin2 * (106) 

This curve represents all the coherent states MQS's of the 
form (|56p . for any value of a and tp. The distance depends 
exclusively from the average number of reflected photons, 
i?|a| 2 , multiplied by a scale factor £ 2 = sin 2 tp. This 
term is proportional to the phase-space distance (FigJSJ) 
between the two components |ae IV )and \ae~ tv ), equal to 
d 2 = 4|a| 2 sin 2 p. In figEHare reported the distances for 
different values of <p for a — 4. For tp = ?, the coherent 
states exhibit opposite phases. Such condition represents 
the limit situation, in which the loss of coherence is higher 
since sin 2 ip = 1. We observe, as is shown by figETJ that 
the value of D is reduced down to ~ 0.1 once i?|a| 2 = 
R(n) = 1. Hence, the loss on the average of a single 
photon cancels most of the coherence in the quantum 
superposition state. 

C. Bures distance for QI-OPA amplified states 

As a following step, we have evaluated numerically 
the distinguishability of {|$ '")} through the distance 
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FIG. 21: (a) Analytical results of the Bures distance for the 
superposition of coherent states, with a = 4 and for differ- 
ent values of the phase, as function of the average number of 
reflected photons, (b) Universal curve of the Bures distance 
for the superposition of coherent states for any a and tp, ob- 
tained by plotting eq. (|106[) as a function of the parameter 
x = R{n)£ v , with f v = sin 2 p. 



D(|$ + ), |$ - )) between the multiphoton states generated 
by QI-OPA. It is found that this property of {|* + ' - )} 
coincides with the MQS Visibility of |\t ), in virtue 
of the phase-covariance of the process: D(\^ + ), \^~)) 
= D(\<$> R ), |$ L ))= D(|$+),|$-». The visibility of the 
MQS {|'3 /+ '~)} has been evaluated numerically analyz- 
ing the Bures distance as a function of the average lost 
photons: x = R < n >. This calculation have been 
performed by taking the complete expression of the den- 
sity matrix, reported in Section IV and Appendices A-B, 
and by performing an approximate calculation of the fi- 
delity through numerical algebraic matrix routines. This 
algorithm has been tested by evaluating numerically the 
Bures distance between the quantum superposition of co- 
herent states | a) ± | — a). The comparison with the ana- 
lytical result of Eq. (|106[) gave a high confidence level for 
the approximate results. The results for different values 
of the gain for equatorial macroqubits are reported in 
FigJ22-(a). 

Note that for small values of x the decay of D{x) 
is far slower than for the coherent state case shown in 
Fig l2U (b). Furthermore, after a common inflection 
point at D ~ 0.6 the slope of all functions D{x) cor- 
responding to different values of < n > increases fast 
towards the infinite value, for increasing x n > 
and: R — * 1. The latter property can be demonstrated 
with a perturbative approach on the density matrix. We 

find, in the low T and high gain limit of — -^^- k - , the 

slope: lim^ lim T _ - Um g _(l + 4C 2 + 

2C 2 r(l + 2r 2 )l?) = co. For total particle loss, R = 1 and 
x = < n > is: D(x) = 0. All this means that the MQS 
Visibility can be significant even if the average number 
x of lost particles is very close to the initial total num- 
ber < n >, i.e. for R ~ 1. This behavior is opposite to 
the case of coherent states where the function D(R\a\ 2 ) 
approaches the zero value with a slope equal to zero: 
FigED(b). 

For sake of completeness, we then performed the same 
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FIG. 22: (a) Numerical evaluation of the distance D(x) 
between two orthogonal equatorial macro- qubits as 
function of the average lost particle x — R < n >. Black 
straight line refers to g = 0.8 and hence to (n) « 4, red dashed 
line to g = 1.1 and (n) « 8, green dash-dotted line to g = 1.3 
and (n) ~ 12. (b) Numerical evaluation of the distance D(x) 
between two orthogonal linear macro-qubits \$ H ' V ) as func- 
tion of the average lost particle x — R < n >. Black straight 
line refers to g = 0.8 and hence to (n) w 4, red dashed line to 
g = 1.1 and (n) « 8, green dash-dotted line to 5 = 1.3 and 
(n) « 12. 



calculation for the linear macroqubits l^^' 17 ), and the re- 
sults are reported in Figl2"2"l-(b). For this injected qubit, 
not lying in the equatorial plane of the Bloch sphere, the 
amplification process does not correspond to an " optimal 
cloning machine" . For this reason, the flow of noise from 
the enviroment in the amplification stage is not the mini- 
mum optimal value, and hence the output states possess a 
faster decoherence rate. Indeed, the output distributions, 
as shown in FigflUl-(b). do not possess the strong unbal- 
ancement in polarization of the equatorial macroqubits 
\<& v ) that is responsible of their resilience structure. 



D. Distinguishability enhancement through an 
Orthogonality-Filter 



As a further investigation, we consider the case of 
a more sophisticated measurement scheme based on an 
electronic device named O-Filter (OF). The demonstra- 
tion of microscopic-macroscopic entanglement by adopt- 
ing the O-Filter based measurement strategy was re- 
ported in [2]]]. The POVM like technique implied 
by this device locally selects, after an intensity measure- 
ment the events for which the difference between the 
photon numbers associated with two orthogonal polar- 
izations — m 7r± I > k, i.e. larger than an adjustable 
threshold k [20] , where is the number of photon polar- 
ized 7r and m 7T± the number of photon polarized 7Tj_ . By 
this method a sharper discrimination between the out- 
put states 1$^) and |$ ¥ ' J -) can be achieved. The action 
of the OF, sketched in Fig[T] can be formalized through 



ma*) = EE", 



mjk) 



n=k m=0 
00 m—k 

m=k n=0 



(107) 

(108) 
(109) 



where the Ii n ^ m = \nn, mir±) (nir, mir± | are Fock-state 
Von-Neumann projectors that describe the performed in- 
tensity measurement. The average of the couple of op- 
erators (F^l(k)) + (H~£l(k)) defines the success prob- 
ability of the O-Filter, i.e. the rate of events leading to 
one of the conclusive outcomes (±1). To calculate the ac- 
tion of the O-Filter in the Bures distance, we projected 
the density matrix of the states over the joint subset cor- 
responding to the (±1) outcomes, neglecting only the 
terms leading to the inconclusive (0) result. Then, the 
same numerical analysis of the previous section has been 
performed. In Fig|23l-fa) the results for g = 0.8 and dif- 
ferent values of k are reported. Note the increase of the 
value of D(x), i.e. of the MQS Visibility, by increasing k. 
Of course, the increase in the MQS Visibility through the 
O-Filter device is achieved at the cost of a lower success 
probability (FigJ23r( a ))- According to the graphical anal- 
ysis of Section HVl on the photon number distributions of 
the equatorial amplified macro-states, the O-Filter de- 
vice improves the MQS Visibility since it exploits the 
peculiar unbalancement in polarization of the equatorial 
amplified macro-states in a Fock-space analysis. In the 
selected regions, the {1$^), |$ ¥ '- L )} states can be discrim- 
inated with a higher fidelity. As a further analysis, we 
reported in figE3r(b) the trend of the success probabil- 
ity as a function of the threshold k for different values 
of the transmittivity. Interestingly enough, we note that 
the success probability depends only on the ratio between 
the threshold k and the number of transmitted photons 
T(n), since the three curves for different T are almost 
superimposed. This is a consequence of the property of 
the photon number distribution of the equatorial macro- 
qubits, since its form is almost left unaltered after the 
propagation over a lossy channel. 



VIII. WIGNER FUNCTION FOR 
NON-DEGENERATE QUANTUM INJECTED 
OPTICAL PARAMETRIC AMPLIFIER IN THE 
EPR CONFIGURATION 

In this section we derive the Wigner function asso- 
ciated to the non-degenerate quantum injected optical 
parametric amplifier, working in a non-collinear EPR 
configuration. The scheme for this device is reported in 

figH 

In the stimulated regime, this amplifier acts as an N — > 
M universal optimal quantum cloning machine (54l . [55| . 
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equations, obtaining: 
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FIG. 23: (a) Numerical evaluation of the Bures distance be- 
tween two orthogonal equatorial macro-qubits after discrimi- 
nation with the O-Filter for different values of the threshold 
k (g — 0.8). The filtering probabilities for the three cases are 
respectively P{k = 4T(n}) = 1.6 * ICT 3 , P{k = GT(n}) = 
7.14 * 10" 5 and P(k = 8T(n}) = 3.06 * 10~ 6 . (b) Success 
probability for the O-Filter for g — 0.8 and different values of 
the transmittivity as a function of the parameter r ^ . The 
latter expresses the threshold k with respect to the average 
number of the incident photons in the detection apparatus. 
We note that the curves for different transmittivities can be 
almost superimposed. 
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In the following paragraphs we shall proceed with the 
calculation of the Wigner function for this amplifier in 
the stimulated regime, i.e. when a generic Fock-state is 
injected on the spatial mode ki. 



A. Wigner function for the non-collinear QIOPA in 
absence of decoherence 

In this section we derive the Wigner function for the 
QIOPA in a non collinear configuration. Without loss 
of generality, let us restrict our attention to the equato- 
rial polarization basis w^, = -j= (tth + e lv nv) and n v>± — 

(— e^^TTH + ttV). As for the collinear case, the in- 
jected state over spatial mode ki is the generic Fock-state 
\ipin) 12 = \N(p, M<p±} 1 (g) \0<p,0<p±)2- The characteristic 
function is then evaluated starting from the definition: 

Xn,m {»?,£,*} = l2 {^ in \e ni& lHW-^H(t) e n24vW-^vW 

^ _(H3) 

Let us apply the transformation {17, — > \T](i), 
where: 

ni = ^= (vi + ^6) ; Ix = (-e^m + 60H4) 



and: 



FIG. 24: Scheme for the non-collinear optical parametric am- 
plification of a single-photon state. An EPR pair is generated 
in a first non-linear crystal. The photon on mode kr acts as a 
trigger to conditionally prepare the photon on mode ki. The 
latter is superimposed with a strong kp pump beam, and is 
amplified in a second non-linear crystal oriented in the same 
optical configuration of the EPR source. The beam-splitters 
on modes ki and k2 in the figure represent the lossy trans- 
mission channel. 



The interaction Hamiltonian of this device can be written 
in the form: 



Hepr = ihx (a\ 7T al 7T± - a\ nx al w j + h.c. (110) 



where n stands for any polarization state, and i = 1,2 
for the two spatial modes. The time evolution of the 
field operators can be directly derived by the Heisenberg 



?7l(2)W = %(2)C- Vl(l)S 
?i(2)M = £r(2)C + ?2(i)^ 



(116) 
(117) 



With analogous calculation to the one performed in Sec- 
tion HIT] we obtain: 



XN,Ai{v,^t} = exp 



4e(i%(*)i 2 +i^(*)i : 



3=1 



N 



^ (N- n)\n\ 2 

\n=0 y ' 



■ M 

E 



M\(-l) m 
(M — m)\m\ 2 



(118) 

The Wigner function of the amplified field can be then 
expressed as the 8-dimensional Fourier transform of the 
characteristic function, according to its definition. By in- 
serting the explicit expression of the characteristic func- 
tion, by separating the integrals on each couple of com- 
plex conjugate variables, and by exploiting the integral 
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relations already used in sec Mil we find: 

W,M^> /?, *} = ) 4 (-1) W+M e - 21 A|2 

L w (|A A {a} + e-^A fl {/3} | 2 ) 
L M (\~A A {a}e^ + A B {/3}| 2 ) 



(119) 



The Aa {a}, Ab {(3} and |A| 2 variables are defined in 
Appendix [Cl where all the details of the calculation of 
the Wigner function are reported. 

Let us we consider the injection of a single photon with 
polarization state \ip) 1 = -A= (\H) + e %v \V)). The Wigner 

function reads [l5j : 

if; 2 

W Mx) P, t} = ^e" 21 A| (| A. 4 {aj+e^As {/?} | 2 - 1) 

(120) 

where |A^ {a} + e '^A^ {/3} | 2 represents the interfer- 
ence term due to the quantum superposition form of 
the input state. Indeed this interference term defines 
a phase-space region where the Wigner function is nega- 
tive, showing the broadcasting of the quantum properties 



J 



of the injected single photon state to the amplified field. 
For more details on the amplification of a single photon, 
refer to fla. 



B. Wigner function for the non-collinear EPR 
optical parametric amplifier in presence of 
decoherence 



In this section we are interested in the Wigner function 
for the non-degenerate optical parametric amplifier in 
presence of losses. The lossy channels on the two spatial 
modes {ki,k2} are again simulated by a beam-splitter 
model, with input-output relations: 



= VTai(i) + iVRb. 



(121) 



where a vacuum state is injected in the 6$ ports of the 
beam-splitters. We considered the case of two transmis- 
sion channels for modes ki and k2 with equal efficiencies. 
The characteristic function is then evaluated as: 



XN,M{r),Z,R,i} = AB ^ in \e^<M-^*' a ^e^ (122) 

The transformation to the {77, £} variables is the same as in the unperturbed case, and has been defined previously in 
Eqq. (|114M115|) . Following the approach of Section IVIII| we insert the beam-splitter relations (I121j) and evaluate the 
averages after writing the exponential operators in normally ordered form. We obtain: 



Xn,m {i],£,R,t} = exp 



~Ul + 2RS 2 )J2(\rj J (t)\ 2 + \H ] (t)\ 
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exp [-RCS fa (t)tj 2 (t) + 1)1 (t)rf 2 (f ))] 
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N\(-T) n 2 „\ / A M\(-T) r ' 
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E 
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(123) 



The Wigner function of the field can then be evaluated as the 8-dimensional Fourier transform: 

{a,(3,R,t} = -^ J J J J f JJ d\d 2 l}j xn.m H e. R, *} exp |e " 

xexp|E[^W^W-^(^W 



\N(p,Mip± 



J = l 



r 



[1241 



After the explicit insertion of the characteristic function, where the {a, R, t} and Ij^ {a, R, t} are: 
the Wigner function can be written as a product of two 
integrals: 



lJf{a,R,t} = -^E d "( T ) J ".» (e',A»iSi(*),5a(t))(126) 

1 M 

lM{P,R,t} = — E<(^)^,o(e',-/x;ft(t)J 2 (a)7) 



Hi 



i} {a, (3, R, t} = 1+ {a, R, t} I M {/?, R, t} 



(125) The J n ,o integrals are explicited in Appendix A in 
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Eq. (|Dl[) . The parameters e d^(T) for the non- 
collincar optical parametric amplifier are: 



N\ (-T) T 
n 



1 

2 

i?C5 



(1 + 2RS 2 



(128) 

(129) 
(130) 



Let us again focus on the single-photon case, corre- 
sponding to N = 1 and M = 0. With an analogous 
procedure to the collinear case, we analyze the persis- 
tence of negative regions in the Wigner function after the 
decoherence process. Analyzing the form of Eq. (|120p . 
we note that the minimum occurs when \Aa {a} + 
e~^A B {/?} | 2 = and |A| 2 = 0. This point corresponds 
to the origin of the 8-dimcnsional phase space given by 
the {a(t), /3(i)} variables. The evolution of the Wigner 
function in this point is explicitly evaluated as: 



U i 



{0,0, R,t} = 



16 



2R-1 



[1 + 4R{1 - R)S 2 ] C 



(131) 



The value in the origin of the Wigner function is re- 
ported in Fig. [25] We note an analogous trend with 
respect to the collinear case, while the absolute value of 
the amount of negativity is smaller. This is due to the 
universal cloning feature of the non-collinear QI-OPA, 
in which the cloning fidelity is smaller than the one in 
phase- covariant case. 

We conclude this Section on the non-collinear QI-OPA 
by stressing that a full phase-space characterization of 
the quantum states generated by this device would re- 
quire a couple of double homodyne measurement setups, 
one for each spatial mode, as the one discussed in Section 

Unci 



IX. CONCLUSION AND PERSPECTIVES 

The quantum properties of the QI-OPA generated 
macroscopic system, have been investigated in phase- 
space by a Wigner quasi-probability function analysis 
when this class of states is transmitted over a lossy chan- 
nel, i.e. in presence of a decohering system - enviro- 
nment interaction which represents the main limitation in 
the implementation of quantum information tasks. We 
first considered the ideal case, in absence of losses, show- 
ing the presence of peculiar quantum properties such as 
squeezing and a non-positive W-representation. Then, 
after a brief review of the properties of the coherent states 
MQS, the resilience to losses of QI-OPA amplified states 
in a lossy configuration was investigated, allowing to ob- 
serve the persistence of the non-positivity of the Wigner 
function in a certain range of the system- enviroment in- 
teraction parameter R. This behaviour was analyzed in 
close comparison with the |a) states MQS, which pos- 
sesses a non-positive W-representation in the same in- 
terval of the interaction parameter R. Moreover, the 
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FIG. 25: Transition from a non-positive distribution to a com- 
pletely positive Wigner function of the single-photon ampli- 
fied states, obtained by the non-collinear QI-OPA in presence 
of decoherence, as a function of the beam-splitter reflectiv- 
ity R. Average photon number of (a) (n) ~ 6.5 and (b) 
(n) ~ 53.5. The negativity is evaluated as the value of the 
Wigner function in the origin. 



more resilient structure of the QI-OPA amplified states 
was emphasized by their slower decoherence rate, repre- 
sented by both the slower decrease in the negative part 
of the Wigner function and by the behavior of the Bu- 
res distance between orthogonal macrostates, the latter 
evaluated in the Fock-Space. Since the negativity of the 
W-representation is a sufficent but not a necessary con- 
dition for the non-classicality of any physical system, fu- 
ture investigations could be aimed to the analysis of the 
decoherence regime in which the Wigner function is com- 
pletely positive, analyzing the presence of quantum prop- 
erties from a different point of view. 

We acknowledge useful discussions with Wolfgang 
Schleich. Work supported by PPJN 2005 of MIUR and 
INNESCO 2006 of CNISM. 



APPENDIX A: DENSITY MATRIX OF THE 
EQUATORIAL AMPLIFIED STATES AFTER 
THE PROPAGATION OVER A LOSSY CHANNEL 



In this appendix we perform the derivation of the den- 
sity matrix of the equatorial amplified states after the 
propagation over a lossy channel. The applied model 
is the same BS-scattering process used throughout the 
whole paper. 

The starting point is the expression of the unperturbed 
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state |$^): 

CO 

l**> = E 



. e ^ — 
C 2 ^ V 2 



6 2 



(Al) 



where C = coshg and T = tanhg, with g gain of the 
amplifier, with the state written in the polarization ba- 
sis {7? v , 7? Vi }. The state can be written in terms of the 
creation operators: 



72 E 



C 2 



rV (4) 2 ' +1 (4J 



(A2) 

We then apply the time evolution independent operators 
Ugg and Ug^\ which describe the interaction in the 
beam-splitter: U BS = e *f ( at6 + a£t ). As a reasonable as- 
sumption, we consider the transmittivity T v and T v± to 
be equal for all polarization states. We can then write 
the evolved state in the Schrodinger picture |4> v ) as: 



\& V )-"BS»BS 



72 X/ 



c 2 



(vt^ + tv^J,) 21+1 (VTcl ± + tVR# v± ) 



-|0) 
(A3) 

We can now exploit the binomial expansion (a + 6)" = 

X)"=o a% b n ~ l {^j > as * ne operators for the two output 

mode of the beam-splitter c, d commute and can be 
treated as algebraic variables. Subsequently applying the 
creation operators to the vacuum we then obtain the fi- 
nal expression of the state after the propagation in the 
beam-splitter: 



i,j=0 

x EE 

m— n— 

n+m 



2i + l 

m 



2i+2j+l-n-m 



x (Vr) (/V/?! 

x |m, n)c ® |2i + 1 — m, 2j — n)o 



(A4) 



The density matrix of the quantum state is then obtained 
as p /lf = |$'^)($'^|. Let us call: 



1 

c 2 



r 

2 



2i + 1 

m 



2i+2j+l-n-m 



Tracing the density matrix //^ on the reflected d BS- 
mode in order to consider our ignorance over the number 
of reflected photons we obtain: 



ft = Tr D (>") = ]T D^.yl*'^) = 
x -y— 

00 00 2»+l 2j 2fe+l 2; 

= E E E E E E^>™"^,p?i m ' n > (p»«i 

® r>(a;,y|2i + 1 - to, 2j - n) D D (2k + l-p,2l- q\x,y) D 

(A6) 

The evaluation of the scalar products between Fock- 
states lead to the desired result. The expression of the 
density matrix on the transmitted c spatial BS-mode is 
then written in the form: 



ft = M^IMIm'MM (A7) 

i,j,k,q=0 



We now report the expressions of the density matrix coefficients, which depend on the parity. For i, j, A:, q even we 



20 



obtain: 
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For i, k odd and j, g even, we obtain: 
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For i, k even and j, g odd, we obtain: 
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Finally, for i,j, k, q odd we obtain: 
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In these expressions, 2-F1 {pi, /?; 7; 2) are Hyper-geometric modes is: 
functions [56| . 



APPENDIX B: DENSITY MATRIX OF THE 

7fff,7fv AMPLIFIED STATES AFTER THE 
PROPAGATION OVER A LOSSY CHANNEL 



The procedure for the evaluation of the density matrix 
of the state after losses is the same applied in the previous 
section. Let us analyze the \& H ) state, whose density 
matrix after the amplification process reads: 



00 n+\ n m+1 m 



c 4 

m,n=0 i=0 j=0 fc=0 9=0 
/ -^\2m+X—k—q 



/ , — \ i+j+k+q / — \ l 

iff) (Vr) 

\iH,jV) c (kH,qV\® 

® \(n + 1 - i)H, (n - j)V) D ((m + 1 - k)H, (m - q)V\ 

(B2) 

After tracing over the reflected mode, we finally obtain: 



n+ 1\ fn\ fm + 1\ fm 
l k ){q 

n-\-m-\-2—i—j — k—q 



n, 771=0 

\(n+l)H,nV) {(m + l)H,mV\ 
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After the insertion of the BS unitary transformation, 
the joint density matrix of the transmitted and reflected 
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where the coefficients "fijk- P are: 
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APPENDIX C: WIGNER FUNCTION FOR THE 
NON COLLINEAR QI-OPA IN ABSENCE OF 
DECOHERENCE 

We begin the calculation of the Wigner function, with- 
out loss of generality, by restricting our attention to the 
equatorial polarization basis, defined by: 
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— e iv tt h + Try) 
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As for the collinear case, the injected state over spa- 
tial mode ki is the generic Fock-state \ip% n ) 12 = 
\Nip, Mip±)i (g) \0ip, 0(p±)2- The characteristic function 
is then evaluated starting from the definition: 

e ^l^ 1V ( t '>-il & iv(t) e f2S,l H {t)~^a,2H{t)^ in ^ i2 

(C3) 

Let us apply the following transformations over the 
{??,£} variables, corresponding to the rotations from the 
{7Tfj,7ry} to the { ifcp , 7?^ j_ } polarization basis: 
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(C5) 
(C6) 
(C7) 



We then explicitly insert the time evolution of the field 
operators (|11HI112|) . Defining: 
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we can re-write the characteristic function in the follow- 
ing form: 



e ?iC t ) a I vi -^C*) s i»>j. e f 2 (t)aI» ; -l2(t)a2 ¥ ,|^. n \ l2 

*" (C12) 

The averages can now be evaluated by writing the op- 
erators in the normally ordered form. With analogous 
calculation to the one performed in Section IIIII we ob- 
tain: 
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(C13) 



The Wigner function of the amplified field can be 
then expressed as the 8-dimensional Fourier transform 
of the characteristic function, according to its defini- 
tion. The integration variables are changed according 



to l , 1 -> {Vj{t),tj(t)} , leading to 
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The calculation then proceeds as follows. By inserting separating the integrals on each couple of complex con- 
the explicit expression of the characteristic function, by jugate variables, and by exploiting the integral relations 
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already used in sec lllll we find: 
W\ Nv , Mvx) {a,P,t} = (1\ (~1) N+M L N (Mai(t)\ 2 ) 

L M (4\A(t)\ 2 ) exp 1-2J2 [l«i(*)| 2 + l^(*)| 2 

(C15) 

As a further step, we can define the following set of 
squeezed and antisqueezed variables: 

1A+ = {a x +a* 2 )e 9 (C16) 

1A- = i{ax-a* 2 )e- 3 (C17) 

7 B+ = (/3i-/3 2 *)e 9 (C18) 

7B _ = z(/3!+/3 2 *)e- s (C19) 

and the following set of quadrature variables: 



f 



a a{«} = -^=(jA+-tJA-) 
As{/3} = (7s+ - «7s-) 



(C20) 
(C21) 



|A| 2 = \ {\lA + \ 2 + \lA-\ 2 + \lB + \ 2 + \lB-[^) 

The Wigner function can then be simply expressed as: 

W ]NvMv±) {a, p, t} = 0) 4 (_i)iv + M e - 2 | 

ijvdA^laj + e-^As^}! 2 ) (C23) 
L M (\-A A {a} e^ + A B {(3}\ 2 ) 

APPENDIX D: INTEGRAL RELATION 

In this appendix we derive the integral relation used 
in the text in Section fVIII Bl We are interested in eval- 
uating integrals of the form: 



J n<m {T,fJ,;z,w) 



2„ j2o i„ \2n\ o\2m 



d ad (3 \a 



=Q * z-az* P'w-Pw* -t(|q| 2 + |/3| 2 ) -^{aP+a* p') 



(Dl) 



We first note that: 

J ntm ( T ,n;z,w) = (-1)"+'™ 

Q2n+2m 



dz n dz* n dw m dw* m 



Jq,o{t,H)Z,w) 



(D2) 



According to this result, it is sufficient to evaluate only 
the Jo,o integral. We start from its explicit expression: 

Jo,o(r,Ai;z,w) = \ [ [ d 2 ad 2 [3e a ' z ~ az * 

71-2 J J (D3) 

e P'w-Pw* e -T(|a| 2 + |/3| 2 ) e -^( a p+a' p') 

We now apply the following | J| = 1 integration variables 
rotation: 



a = cos 8 7 + sin 8 S 
(3 = — sin 8 7 + cos 8 S 



(D4) 
(D5) 



Choosing cos 8 = sin 8 = ;4= we obtain: 
J ,o(r,wz,«;) = i J J d^dHe^-^ 

e rw-5«J* e -r(| 7 | 2 + |5| 2 ) e -i M (5 2 +5* 2 - 7 2 - 7 - 2 ) 

where we defined the rotated parameters: 

1 



(D6) 



V2 
1 



(z — w) 



w = — {z + w) 



(D7) 
(D8) 



The integrals over 7 and 5 are now uncoupled, and 
can be evalated separately by using Eq. (|75|) . Exploiting 
this integral relation, and applying the inverse rotation 
{z, w} — > {z,w}, we obtain: 



Jo,o(t,(j,;z,w) = 



1 



~2 2 eXP 

t — 



2 2 



exp 



T 2 — /i 2 



J (D9) 



The latter result allow then to explicitly evaluate the 
Jn,m integrals, by starting from Jo,o an d performing the 
opportune derivatives, according to Eq. (|D2[) . 
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